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ABSTRACT 

Haney-type  surface  thermal  boundary  conditions  linearly  connect  net  downward  surface  heat  flux  Q  to  air- 
sea  temperature  difference  (gradient-type  condition)  AT,  or  to  climate/synoptic  sea  temperature  difference  (re- 
storing-type  condition)  A T2  by  a  coupling  coefficient  k.  In  this  study,  the  authors  used  the  global  reanalyzed 
data  (6-h  resolution)  of  Q,  surface  air  temperature  TA,  and  sea  surface  temperature  T0  from  the  National  Centers 
for  Environmental  Prediction  during  1  October  1994-31  December  1995  to  verify  the  validity  of  Haney-type 
surface  thermal  boundary  conditions.  First,  daily  means  of  these  variables  were  computed  to  get  rid  of  diurnal 
variation.  Second,  the  cross-correlation  coefficients  (CCC)  between  Q  and  (Arl5  A7\)  were  calculated.  The 
ensemble  mean  CCC  fields  show  (i)  no  correlation  between  Q  and  A T2  anywhere  in  the  world  oceans,  (ii)  no 
correlation  between  Q  and  AT,  in  the  equatorial  regions,  and  (c)  evident  correlation  (CCC  >  0.7)  between  Q 
and  AT,  in  the  middle  and  high  latitudes.  Third,  the  variance  analysis  was  conducted  and  a  value  of  70  W  m~2 
K  1  (65  W  m-2  K  _1)  was  suggested  for  the  coupling  coefficient  k  in  the  northern  (southern)  middle  and  high 
latitude  zone. 

Thus,  the  authors  find  that  the  restoring-type  surface  thermal  conditions  by  no  means  represent  the  net  air- 
ocean  heat  flux  anywhere  in  the  world  oceans.  However,  the  gradient-type  surface  thermal  condition  represents 
the  net  heat  flux  quite  well  for  the  middle  and  high  latitudes.  In  addition,  it  is  also  found  that,  if  the  solar 
shortwave  component  is  treated  separately,  the  gradient-type  condition  will  have  more  fidelity  for  the  middle 
and  high  latitudes. 


1.  Introduction 

The  ocean  is  thermally  driven  by  the  net  downward 
flux  of  heat  across  the  ocean  surface,  Q,  which  is  the 
sum  of  the  downward  flux  of  solar  radiation,  Rs,  minus 
the  net  upward  flux  of  longwave  (or  “back”  radiation) 
Rb,  sensible  heat  QH,  and  latent  heat  QE, 

Q  =  Rs  -  (Rb  +  Q„  +  Qe).  (1) 

Under  the  assumption  that  the  ocean  is  in  contact  with 
an  atmospheric  equilibrium  state  (i.e.,  an  atmosphere 
with  a  near-infinite  heat  capacity),  Haney  (1971)  ob¬ 
tained  a  very  simple  heat  flux  formulation, 

Q  =  kAT,  AT  =  T*  -  T0,  (2) 

where  T%  is  the  apparent  atmospheric  equilibrium  tem¬ 
perature,  Ta  the  sea  surface  temperature  (SST),  and  k 
is  a  coupling  coefficient.  In  spite  of  its  temporal  vari¬ 
ation,  the  parameter  is  usually  taken  as  positive  values 
between  10  and  50  W  m~2  Kr1  by  various  authors  (e.g., 
Marotzke  1994;  Cai  and  Chu  1996).  Here  T*  should  be 
computed  from  surface  heat  fluxes  and  their  dependence 


Corresponding  author  address:  Prof.  Peter  C.  Chu,  Department  of 
Oceanography,  Naval  Postgraduate  School,  Monterey,  CA  93943- 
5000. 

E-mail:  chu@nps.navy.mil 


on  temperature,  as  done  by  Han  (1984)  and  Oberhuber 
(1988).  Despite  Haney’s  emphasis  that  Q  is  proportional 
to  7’*  —  Ta,  many  subsequent  investigators  have  re¬ 
placed  7’*  by  Ta  or  T'%.  Here  TA  is  the  surface  air  tem¬ 
perature  (SAT),  and  T%  is  either  the  climatological  mean 
( T0 )  or  observed  SST.  In  this  paper,  we  verify  the  case 
where  T%  =  (T0).  Thus,  the  Haney  condition  becomes 

Q  =  kA Tu  AT \  =  Ta-  T0,  (3) 

which  is  the  gradient-type  condition,  and 

Q  =  T2,  AT2  =  (Ta)  -  T0,  (4) 

which  is  the  restoring-type  condition.  Pierce  (1996)  point¬ 
ed  out  two  distinct  disadvantages  of  the  restoring-type 
conditions:  (i)  it  imposes  a  time  lag  in  the  model’s  surface 
tracer  field  and  (ii)  it  systematically  underrepresents  the 
model’s  surface  variability.  Under  the  assumption  that  the 
surface  fluxes  contribute  more  to  SST  changes  than  do 
internal  ocean  processes.  Pierce  proposed  a  method  to  cor¬ 
rect  the  amplitude  and  phase  of  T0  for  (4). 

Another  trend  is  to  separate  the  solar  radiation  from 
the  net  downward  flux  (e.g.,  Blumberg  and  Mellor  1987) 

Q  =  Rs  +  Q,  Q=  -(R„  +  Qh+  Qe)  (5) 

and  to  replace  Q  by  Q  in  the  gradient-type  surface 
boundary  condition 


Q  =  kAT1 


(6) 
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and  the  restoring-type  surface  boundary  condition 

Q  =  kA  T2.  (7) 

We  call  Q  the  net  downward  heat  flux  with  solar  radi¬ 
ation,  and  Q  the  downward  heat  flux  without  solar  ra¬ 
diation. 

Since  the  Haney-type  boundary  conditions  have  been 
used  by  the  community  for  more  than  two  decades,  it 
is  time  to  verify  the  validity  of  these  conditions.  How 
good  are  these  conditions?  Is  there  any  temporal  and 
spatial  variation  of  the  validity  of  these  conditions? 

2.  Statistical  evaluation 

a.  Dataset 

We  chose  the  National  Centers  for  Environmental 
Prediction  (NCEP)  reanalyzed  global  surface  fluxes  Rs, 
Rh,  Qh,  Qe,  and  temperatures  Ta  (skin  temperature),  T A 


data  on  a  1.875°  X  1.875°  grid  for  this  study.  The  tem¬ 
poral  resolution  of  the  original  data  is  6  h.  Since  Haney- 
type  surface  thermal  boundary  conditions  are  commonly 
used  for  ocean  modeling  on  a  seasonal  timescale  or 
longer,  we  first  computed  daily  means  of  these  variables 
to  filter  out  diurnal  variation  and  to  obtain  a  dataset  of 
global  Q(x,  y,  d),  Q(x,  y,  d),  TA(x,  y,  d),  Ta(x,  y,  d)  for 
each  day,  d,  from  1  October  1994  to  31  December  1995. 
Then,  from  Ta(x,  y,  d)  data,  we  computed  the  monthly 
mean  values  (T0)  at  each  grid. 

b.  Cross-correlation  coefficient 

The  validity  of  Haney-type  boundary  conditions  (3), 
(4),  (6),  and  (7)  should  be  first  tested  by  cross-corre¬ 
lation  coefficients  (CCC)  between  the  two  time  series: 
Q  [ Q  or  Q\  and  A T  [ATj  or  AT2]  at  day  d,  which  is 
computed  in  this  paper  from  a  91 -day  subset  between 
45  days  prior  to  and  45  days  after  that  date, 


S  {[Q(x,  y,  t,)  -  Q(x,  y)][A7Xx,  y,  r,)  -  AT(x,  y)]} 

Rq,aM)  =  — - — - 7 7 -  (d  -  45  <  T,.  <  d  +  45),  (8) 


where  N  =  91  days,  (x,  y)  the  horizontal  coordinates, 
Q{x,  y)  and  A7’(x,  y)  the  temporal  means  of  the  subset 
data, 

—  1  N 

Q(x,  y)  =  ~  2  Q(x,  y,  t), 

TV  ,=i 

A T(x,  y)  =]-Y  AT(x,  y,  r,),  (9) 

i=  i 

and  crQ(d)  and  crAT(d)  the  standard  deviations  of  the 
subset  data. 


c.  Determination  of  k  by  no-missing  variance 

Since  the  value  of  k  is  quite  uncertain,  such  a  variance 
comparison  between  ( Q ,  Q )  and  (ATt,  A T2)  depends  on 
the  choice  of  k.  On  the  other  hand,  we  may  determine 
the  value  of  k  by  no-missing  variance 

a0{d) 

K(d)  =  -^TfRasT  (ID 

aAT(d) 

and  check  if  the  values  are  within  a  reasonable  range. 

3.  Statistical  tests 


fXQ{d)  =  Jff— J  S  \.Q(x,  y,  t)  ~Q{x,  y)]2 


crAT(d)  =  2  [A7'(x,  y,  t;)  -  A T(x,  y)]2 .  (10) 

Thus,  we  established  little  over  one  year  (368  days: 
16  November  1994-18  November  1995)  four  CCC  da¬ 
tasets:  RaATl,  Rq,ATi ,  Rq.stp  and  RaATl.  The  higher  the 
values  of  CCC,  the  higher  the  confidence  of  a  linear 
relationship  between  ( Q ,  Q)  and  (A7\,  A T2). 

Besides  the  correlation  analysis,  we  still  need  to  do 
the  variance  analysis,  because  it  does  matter  if  ( Q ,  Q) 
and  (A7\,  A7’2)  correlate  closely,  but  typically  used  val¬ 
ues  of  k  (10-50  W  m~2  Kr1)  end  up  missing  all  the 
variance. 


Whether  the  sample  values  of  RQ  AT  represent  good 
or  bad  linear  relationships  (3),  (4),  (6),  and  (7)  between 
(G>  Q)  and  (AT1,,  A Tf)  should  be  tested.  First,  we  use 
the  1  test  to  identify  the  region  with  small  values  of  CCC 
in  which  there  is  no  correlation  between  (Q,  Q )  and 
(AT,,  A r2).  Second,  we  set  up  a  prior  significant  value 
(pu)  for  CCC,  and  use  the  z  test  to  see  where  the  values 
of  Rq  at  are  significantly  larger  than  p0. 


a.  No  correlation  between  Q  and  AT  (t  test) 


The  t  value  for  any  location  (x,  y),  computed  by 


Rb,atVN  -  2 

Vi  -  rqa, 


(12) 


is  used  to  assess  the  statistical  significance  of  RQ  AT  dif- 
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Gradient-type  with  solar  radiation  (ccc) 


Gradient-type  without  solar  radiation  (ccc) 


Fig.  1.  Ensemble  mean  cross-correlation  coefficients  between  A T1  and  net  surface  heat  flux:  (a)  with  solar 
radiation  Q,  (b)  without  solar  radiation  Q,  and  (c)  correlation  coefficient  (X100)  between  monthly  anomalous 
latent  and  sensible  heat  flux  and  changes  in  SST  anomalies  (Fig.  lc  from  Cayan  1992). 


ferent  from  0.  We  begin  with  the  null  hypothesis  that 
the  two  variables  Q  and  AT  are  not  positively  correlated. 
The  significance  level  a  is  the  probability  that  the  given 
value  of  1  is  exceeded  purely  by  chance.  Equivalently, 
it  is  the  probability  of  incorrectly  rejecting  the  null  hy¬ 
pothesis.  For  a  given  significance  level  a  (a  =  0.05 
was  used),  we  compare  the  t  value  with  ta:  if  t  <  ta, 
we  accept  the  null  hypothesis,  and  there  is  no  linear 
(positive)  relationship  between  Q  and  AT  at  the  given 
location  ( x ,  y). 


b.  Significant  correlation  (RQAT  >  p0)  between  Q  and 
AT  (z  test ) 


The  z  value,  computed  by 

(1  +  ^o,Ar)(l  ~  Po) 

(1  -  *flAr)(l  +  Po) 
is  used  to  assess  the  statistical  significance  of 

Rq.at  Po 


Vn~ 


z  = 


In 


(13) 


(14) 
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Fig.  1.  (Continued) 
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for  any  location  (x,  y) .  This  z  value  satisfies  the  normal 
distribution. 

We  start  with  the  null  hypothesis  that  H0 :  RQAT  =  p0, 
H1:  Rq  at  >  p0.  The  significance  level  a  is  the  proba¬ 
bility  that  the  given  value  of  z  is  exceeded  purely  by 
chance.  For  a  given  significance  level  a  {a  =  0.05  was 
used),  we  compare  the  z  value  with  za :  if  z  >  za,  we 
accept  H1,  and  the  CCC  between  Q  and  AT  exceeds  p0 
at  the  given  location  (x,  y).  In  this  study,  we  use  p0  = 
0.7. 

After  RqATi(x,  y,  d),  RaATt (x,  y,  d ),  RaATl(x,  y,  d), 
and  RqATi(x,  y,  d)  are  calculated,  we  used  (12)  to  com¬ 
pute  t  values  and  (13)  to  calculate  the  z  values  for  each 
grid  point,  compared  the  t  values  to  ta  to  test  if  (Q,  Q) 
and  (AT1!,  A T2)  were  independent,  and  compared  the  z 
values  to  za  to  test  if  RQA r,(x,  y,  d),  RqA r,(x,  y,  d), 
Rq.at2(x,  y,  d),  and  Rq.atXx,  y,  d)  were  significant. 

4.  CCC  fields  over  the  world  oceans 

Four  time  series,  RQA Ti(x,  y,  d),  Rq:ATi{x,  y,  d), 
RQAt2(x,  y,  d),  and  Rq.atXx,  y,  d),  were  decomposed 
into  ensemble  means: 

<*e>An>  —  2  Rq.atM’  I7’ 

(Rq.at,)  =  2  Rq,ai\ (x>  y>  d), 

( Rq.AT2 )  —  TTn  2  Rq,AT2(X>  y>  d), 

JOo  d 

(Rq.at2)  =  77T  2  Rq.atSx,  y,  d),  (15) 

30o  d 

which  show  the  overall  linear  relationships  between  (Q, 
Q )  and  (A 7) ,  AT,)  and  anomalies: 


Rq.at1  —  Rq,at2  ~~  (Rq,at1)-’ 

Rq,at1  —  Rq,at2  (Rq.at'I’ 

Rq.at2  ~~  Rq,at2  ( Rq.at 

R’q.at2  =  Rq,at2  ~  (Rq.at)’  (16) 

which  indicate  temporal  variations  of  such  relationships. 

a.  Ensemble  mean  CCC  fields 

1)  Gradient- type  conditions 

Figures  la  and  lb  illustrate  the  global  distributions 
of  {Rq.at)  and  < RqAT ,),  respectively.  The  region  with 
large  values  (sO.7)  of  (RQAT)  (or  ( RqAT ,))  is  called  the 
valid  region,  and  the  region  with  small  values  (^0.2) 
of  ( Rq.at ,)  (or  ( Rq,at ,))  is  called  the  invalid  region.  The 
gradient-type  Haney  boundary  conditions  with  and 
without  solar  radiation  are  valid  in  middle  and  high 
latitudes  (higher  than  30°)  away  from  coasts,  and  are 
invalid  in  equatorial  regions  (10°N-10°S).  The  invalid¬ 
ity  in  the  equatorial  and  coastal  regions  may  imply  the 
improper  representation  of  latent  heat  QE  by  temperature 
difference  A Tl. 

Under  the  assumption  that  the  upper  ocean  is  a  con¬ 
stant-depth  well-mixed  slab  that  exchanges  heat  with 
the  atmosphere,  Cayan  (1992)  computed  CCC  (Fig.  lc) 
between  monthly  anomalous  latent  and  sensible  heat 
flux  and  changes  in  SST  anomalies  over  the  North  At¬ 
lantic  and  North  Pacific  from  the  monthly  summaries 
trimmed  subset  of  the  Comprehensive  Ocean-Atmo¬ 
sphere  Data  Set  (COADS).  Comparison  of  Fig.  la  and 
lb  with  Fig.  lc  shows  a  substantial  agreement  (high 
CCC  in  middle  and  high  latitudes  and  low  CCC  in  equa¬ 
torial  regions). 
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Restoring-type  with  solar  radiation  (ccc) 
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Restoring-type  without  solar  radiation  (ccc) 


Fig.  2.  Ensemble  mean  cross-correlation  coefficients  between  A7\  and  net  surface  heat  flux  (a)  with  solar 
radiation  Q  and  (b)  without  solar  radiation  Q. 


2)  Restoring-type  CONDITIONS 

Figure  2  illustrates  the  global  distributions  of 
(Rq  ATi)  and  (RqA12),  respectively.  The  small  values  (less 
than  0.2)  of  ( RQAT ,)  and  ( RqATi )  are  all  over  the  globe 
except  the  eastern  equatorial  Pacific  and  South  Atlantic 
section  near  Antarctica,  where  (Rq  Ar,)  has  values  around 
0.3.  No  correlation  between  ( Q ,  Q )  and  A T2  may  in¬ 
dicate  that  the  restoring-type  surface  thermal  condition 
by  no  means  represents  the  net  air-ocean  heat  flux  any¬ 
where  in  the  world  oceans. 


b.  Anomaly  CCC  fields 

We  used  empirical  orthogonal  function  (EOF)  anal¬ 
ysis  to  delineate  the  major  modes  of  variability  of  the 
anomaly  CCC  fields:  R'a  Ar,  and  R'a ATi,  which  are  de¬ 
composed  by 

R'  =  2  p CJ.d)K(x»  yj),  (17) 

a 

where  PC  a(d)  is  the  principal  component,  representing 
the  temporal  variation  of  the  associated  spatial  pattern 
described  by  EOF  </>„(x,,  y;).  The  time  series  is  anal- 
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Table  1.  Variances  of  the  first  six  leading  EOFs  for  the  gradient- 
type  conditions. 


EOF 

With  solar 
radiation 

Without  solar 
radiation 

1 

0.643 

0.742 

2 

0.179 

0.158 

3 

0.073 

0.040 

4 

0.038 

0.020 

5 

0.017 

0.009 

6 

0.012 

0.009 

ogous  to  a  projection  of  CCC  anomaly  through  the  “fil¬ 
ter”  of  an  EOF  mode  during  a  timescale. 

The  first  two  leading  EOFs  are  able  to  account  for 
82%  (for  and  90%  (iovR'g  ATl)  of  the  total  variance 
during  the  period  (Table  1).  Each  EOF  mode  is  nor¬ 
malized  so  that  its  spatial  variance  is  equal  to  unity.  So, 
the  patterns  of  the  first  two  EOFs  are  enough  to  explain 
the  spatial  anomalies  of  the  global  CCCs.  Here,  we  show 
only  EOF1  and  PC,  for  R  as  an  example. 

The  EOF1  mode  (Fig.  3a)  has  a  dipole  pattern  fea¬ 


turing  the  Northern  Hemisphere  positive  (maximum  val¬ 
ue  near  0.02)  and  Southern  Hemisphere  negative  (min¬ 
imum  value  near  —0.02),  and  accounts  for  up  to  64.3% 
of  the  total  spatial  variance.  This  pattern  is  believed  to 
be  related  to  the  solar  radiation.  The  first  principal  com¬ 
ponent,  PC, (d),  during  the  integration  period  is  shown 
in  Fig.  3b.  EOF1  mode  <£,(xf,  y.)  is  generally  positive 
(negative)  in  the  Northern  (Southern)  Hemisphere  (Fig. 
3a).  Therefore,  PC,(g()  >  0  corresponds  to  positive  (neg¬ 
ative)  CCC  anomalies  in  the  Northern  (Southern)  Hemi¬ 
sphere,  and  PC,(af)  <  0  corresponds  to  negative  (pos¬ 
itive)  CCC  anomalies  in  the  Northern  (Southern)  Hemi¬ 
sphere.  Figure  3b  shows  a  strong  seasonal  variation: 
From  late  October,  when  PC,  (c/)  =  0,  it  increases  with 
time  until  mid- January  when  PC,  reaches  the  maximum 
value  of  9.0,  and  then  decreases  to  0  in  early  March. 
After  early  March,  PC,  becomes  negative  until  late  Oc¬ 
tober  and  reaches  the  minimum  value  of  —7.0  in  early 
July.  From  late  October  to  early  March,  the  Northern 
(Southern)  Hemisphere  is  winter  (summer);  PC, (of)  > 
0.  The  contribution  of  EOF1  to  the  total  CCC  field. 


EOF  1  Gradient-type  with  solar  radiation 


Restoring-type  with  solar  radiation  PC  1 


Fig.  3.  (a)  First  EOF  mode  of  R'QATl  and  (b)  associated  time  series  of  PC,. 
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Gradient-type  with  solar  radiation 


month  /  date 


Gradient-type  without  solar  radiation 


month  /  date 


Fig.  4.  Seasonal  variation  of  cross-correlation  coefficients  averaged 
within  five  geographic  zones  between  (a)  Q  and  AY,  and  (b)  Q  and 
A F, .  Here  the  thin  solid  curve  indicates  the  southern  middle  and  high 
latitude  zone  (south  to  30°S),  the  thick  solid  curve  indicates  the  north¬ 
ern  middle  and  high  latitude  zone  (north  to  30°N),  the  dotted  curve 
indicates  the  southern  subtropical  zone  (10°-30°S),  the  dashed-dotted 
curve  indicates  the  northern  subtropical  zone  (10°-30°N),  and  the 
dashed  curve  indicates  the  equatorial  zone  (10°S-10°N). 

PC, (of)  X  </>,(x,,  yj),  adds  positive  values  to  the  ensem¬ 
ble  CCCs  with  a  maximum  increasing  value  of  0.18 
(  =  9.0  X  0.02)  in  the  Northern  Hemisphere  and  adds 
negative  values  to  the  ensemble  CCCs  with  a  maximum 
decreasing  value  of  —0.18  [=9.0  X  (—0.02)]  in  the 
Southern  Hemisphere.  On  the  other  hand,  from  early 
March  to  late  October  PC,(g0  <  0.  The  contribution  of 
EOF1  to  the  total  CCC  field  is  reversed  to  the  period 
from  late  October  to  early  March.  Thus,  the  winter  hemi¬ 
sphere  has  better  cross  correlations  than  the  summer 
hemisphere.  The  maximum  temporal  variability  can 
reach  0.36  (=2  X  0.18).  This  is  because  the  surface 
solar  radiation  cannot  be  well  represented  by  the  tem¬ 
perature  difference  AT,. 

5.  Seasonal  CCC  variation  in  five  geographic  zones 

On  the  basis  of  the  spatial  variation  of  the  ensemble 
mean  field,  we  divided  the  world  oceans  into  five  geo¬ 


Restoring-type  with  solar  radiation 


month  /  date 


Restoring-type  without  solar  radiation 


Fig.  5.  Seasonal  variation  of  cross-correlation  coefficients  averaged 
within  five  geographic  zones  between  (a)  Q  and  A T2  and  (b)  Q  and 
A T2.  Here  the  thin  solid  curve  indicates  the  southern  middle  and  high 
latitude  zone  (south  to  30°S),  the  thick  solid  curve  indicates  the  north¬ 
ern  middle  and  high  latitude  zone  (north  to  30°N),  the  dotted  curve 
indicates  the  southern  subtropical  zone  (10°-30°S),  the  dashed-dotted 
curve  indicates  the  northern  subtropical  zone  (10°-30°N),  and  the 
dashed  curve  indicates  the  equatorial  zone  (10°S-10°N). 


graphic  zones:  equatorial  zone  (10°N-10°S),  northern 
subtropical  zone  (10°-30°N),  southern  subtropical  zone 
(10°-30°S),  northern  middle  and  high  latitude  zone 
(north  of  30°N),  and  southern  middle  and  high  latitude 
zone  (south  of  30°S).  At  each  day,  we  averaged  the  CCC 
data  RaATi(x,  y,  d ),  Rq,ATi(x,  y,  d),  RQATl(x,  y,  d),  and 
RqATXx<  d)  spatially  within  each  zone.  Thus,  we  ob¬ 
tained  five  time  series  for  each  of  the  CCC  data: 

Rh,AT,(d),  RsQAT2(d),  RsaATl(d).  Here,  the  su¬ 
perscript  “5”’  means  the  spatial  average. 

a.  Gradient-type  conditions 
1)  With  solar  radiation 

Figure  4a  illustrates  the  seasonal  variation  of  RsgATi(d) 
for  the  five  geographic  zones.  The  value  fluctuates  from 
—0.1  to  0.05  in  the  equatorial  zone  (no  correlation  be- 
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Gradient-type  with  solar  radiation 


Gradient-type  without  solar  radiation 


Fig.  6.  Seasonal  variation  of  coupling  coefficient  k  averaged  within 
five  geographic  zones  determined  by  no  loss  of  variance  for  (a)  Q 
=  kATx  and  (b)  Q  =  kATx.  Here  the  thin  solid  curve  indicates  the 
southern  middle  and  high  latitude  zone  (south  to  30°S),  the  thick  solid 
curve  indicates  the  northern  middle  and  high  latitude  zone  (north  to 
30°N),  the  dotted  curve  indicates  the  southern  subtropical  zone  (10°— 
30°S),  the  dashed-dotted  curve  indicates  the  northern  subtropical 
zone  (10°-30°N),  and  the  dashed  curve  indicates  the  equatorial  zone 
(10°S-10°N). 


tween  Q  and  A Tt),  from  0.22  to  0.54  in  the  southern 
subtropical  zone,  and  from  0.1  to  0.54  in  the  northern 
subtropical  zone  (weak  correlation  between  Q  and  A 7’, ). 
However,  the  value  fluctuates  from  0.7  to  0.86  all  year 
round  in  the  southern  middle  and  high  latitude  zone 
(high  correlation),  and  from  0.5  to  0.7  in  the  summer 
half-year  (23  March-1  October)  and  from  0.7  to  0.86 
in  the  winter  half-year  (1  October-23  March)  in  the 
northern  middle  and  high  latitude  zone.  On  the  basis  of 
correlation,  the  surface  thermal  boundary  condition,  Q 
=  kAT{,  is  valid  only  in  the  middle  and  high  latitudes. 
The  cross-correlation  coefficient  RsQATl{d)  has  seasonal 
variation  in  both  northern  and  southern  middle  and  high 
latitude  zones:  higher  values  during  winter  and  lower 
values  during  summer.  This  may  be  caused  by  improper 
representation  of  the  solar  radiation  Rs  by  the  temper¬ 
ature  difference  AT). 


Restoring-type  with  solar  radiation 


Restoring-type  without  solar  radiation 


month  /  date 


Fig.  7.  Seasonal  variation  of  coupling  coefficient  k  averaged  within 
five  geographic  zones  determined  by  no  loss  of  variance  for  (a)  Q 
=  kAT2  and  (b)  Q  =  kA7\.  Here  the  thin  solid  curve  indicates  the 
southern  middle  and  high  latitude  zone  (south  to  30°S),  the  thick  solid 
curve  indicates  the  northern  middle  and  high  latitude  zone  (north  to 
30°N),  the  dotted  curve  indicates  the  southern  subtropical  zone  (10°— 
30°S),  the  dashed-dotted  curve  indicates  the  northern  subtropical 
zone  (10°-30°N),  and  the  dashed  curve  indicates  the  equatorial  zone 
(10°S~10°N). 

2)  Without  solar  radiation 

Figure  4b  illustrates  the  seasonal  variation  of7?|Ar  (d) 
for  the  five  geographic  zones.  The  correlation  is  higher 
between  Q  and  A 7)  than  between  Q  and  A 7),  indicating 
the  benefit  of  excluding  the  solar  radiation  from  the  net 
surface  heat  flux.  For  example,  RSQATl(d)  is  near  0.86 
year  round  in  the  southern  middle  and  high  latitudes, 
and  varies  between  0.74  and  0.84  in  the  northern  middle 
and  high  latitudes.  However,  R^AI)(d)  is  still  very  low 
(between  0  and  0.2)  in  the  equatorial  zone  and  becomes 
higher  but  less  than  0.7  in  both  northern  and  southern 
subtropical  zones. 

b.  Restoring-type  conditions 

Figure  5  illustrates  the  seasonal  variations  of7?^Ar]  (d) 
and  RSQATl(d)  for  the  five  geographic  zones,  respectively. 
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Gradient-type  without  solar  radiation  (ccc) 


Gradient-type  without  solar  radiation  (ccc) 


Fig.  8.  Ensemble  mean  cross-correlation  coefficients  between  AY,  and  net  surface  heat  flux  without  solar 
radiation  Q  after  the  daily  NCEP  dataset  undergoes  (a)  30-day  running  mean,  (b)  60-day  running  mean,  and 
(c)  90-day  running  mean. 


No  matter  where  and  when  it  is,  there  is  no  correlation 
between  ( Q ,  Q )  and  A T2  in  the  world  oceans. 

6.  Seasonal  variation  of  coupling  coefficient  k 

We  computed  the  coupling  coefficient  k  using  (11) 
and  then  averaged  k  spatially  within  each  of  the  five 
zones.  Thus,  we  obtained  both  the  spatial  (five  zones) 
and  the  temporal  variability  of  k.  Theoretically,  k  can 
be  time  dependent.  Practically,  k  is  usually  taken  as  a 
constant  in  the  surface  thermal  boundary  conditions  (3), 
(4),  (6),  and  (7).  Therefore,  besides  high  cross  corre¬ 


lations,  we  should  also  consider  quasi-steadiness  of  k 
in  evaluating  the  surface  thermal  boundary  conditions. 

a.  Gradient  type 

For  surface  net  heat  flux  with  solar  radiation  (Fig. 
6a),  k  varies  from  —10  to  10  W  m~2  Ku1  in  the  equa¬ 
torial  zone,  from  15  to  78  W  m  2  K  1  in  the  southern 
subtropical  zone,  from  36  to  80  W  m  2  K  1  in  the  north¬ 
ern  subtropical  zone,  from  70  to  95  W  m  2  K  1  in  the 
southern  middle  and  high  latitude  zone,  and  from  60  to 
70  W  m  2  K  1  in  the  northern  middle  and  high  latitude 


May  1998 


CHU  ET  AL  . 


899 


Gradient-type  without  solar  radiation  (ccc) 


Longitude  (E),  Dep?h=0(m) 

Fig.  8.  ( Continued ) 


zone.  When  the  solar  radiation  is  excluded  from  the  net 
surface  heat  flux  (Fig.  6b),  k  has  large  temporal  vari¬ 
ation  in  the  equatorial  zone.  However,  k  has  small  tem¬ 
poral  variation  (56-74  W  m~2  Kr1)  in  the  middle  and 
high  latitudes  of  both  Northern  and  Southern  Hemi¬ 
spheres. 

On  the  basis  of  (i)  high  correlation  between  Q  and 
AT, ,  and  (ii)  quasi-steadiness  of  k,  we  may  conclude 
that  the  surface  thermal  boundary  condition 

Q  =  «AT , 

is  a  good  parameterization  of  surface  net  heat  flux  for 
the  middle  and  high  latitudes. 

b.  Restoring  type 

The  coupling  coefficient  k  has  large  spatial  and  tem¬ 
poral  variation  (from  —70  to  25  W  m  2  K  1 )  for  both 
including  and  excluding  solar  radiation  (Fig.  7).  Several 
studies  (Marotzke  1994;  Rahmstorf  and  Willibrand 
1995;  Seager  et  al.  1995;  Pierce  et  al.  1996;  Marotzke 
and  Pierce  1997)  show  the  lengthscale  dependence  of 
values  of  k  for  the  restoring-type  boundary  condition. 
For  example,  Seager  et  al.  (1995)  found  that  for  typical 
anomaly  length  scales,  a  value  of  k  of  15  W  m~2  K~' 
is  appropriate. 

Since  (i)  low  correlation  between  ( Q ,  Q)  and  A T2, 
and  (ii)  rapid  variation  of  k,  we  may  conclude  that  the 
restoring-type  conditions  do  not  represent  any  surface 
thermal  forcing  and  in  turn  any  discussion  on  k  for  this 
type  of  condition  is  useless. 

7.  Timescale  dependence  of  k  for  the  gradient-type 

conditions 

The  coupling  coefficient  k  has  units  of  watts  per 
square  meter  per  kelvin.  A  timescale  can  be  derived  as 


t  =  pcph/K,  where  p  is  the  density,  cp  is  the  specific 
heat  of  seawater,  and  h  is  a  typical  mixed  layer  depth. 
Weaver  and  Sarachik  (1991)  used  a  timescale  of  25  days. 
Marotzke  and  Willibrand  (1991)  used  30  days.  For  a 
mixed  layer  50  m  deep  these  correspond  to  values  of  k 
of  97  W  m~2  K  1  and  81  W  m~2  K  ',  respectively.  As 
pointed  out  by  a  number  or  authors  (Zhang  et  al.  1993; 
Power  et  al.  1994;  San  ter  et  al.  1995;  Rahmstorf  and 
Willibrand  1995;  Pierce  et  al.  1995;  Cai  and  Chu  1996), 
change  of  values  of  k  has  detrimental  effects  on  the 
realism  of  a  modeled  thermohaline  circulation.  Thus, 
use  of  a  proper  k  value  becomes  an  important  issue. 
Section  6a  shows  that  k  varies  from  60  to  76  W  m~2 
K  1  in  the  middle  and  high  latitudes  of  both  Northern 
and  Southern  Hemispheres  based  on  the  daily  data.  The 
gradient-type  boundary  conditions  are  used  in  ocean 
models  when  the  variation  is  on  different  timescales. 
Many  of  them  are  on  a  seasonal  timescale.  In  order  to 
investigate  the  timescale  dependence  of  k  values,  we 
use  30-day,  60-day,  and  90-day  running-mean  operators 
on  the  original  NCEP  Rs,  Rb,  QH,  QE,  T0,  and  TA  to 
filter  out  high-frequency  variations  and  obtain  three  da¬ 
tasets  of  global  Q{x,  y,  d ),  Q(x,  y,  d),  TA(x,  y,  d),  and 
T0(x,  y,  d),  respectively.  The  smoothed  datasets  have 
different  durations:  16  October  1994-16  December 
1995  for  a  30-day  smoothing,  31  October  1994-1  De¬ 
cember  1995  for  a  60-day  smoothing,  and  15  November 
1994-16  November  1995  for  a  90-day  smoothing.  For 
each  smoothed  dataset,  we  use  (8)  and  (11)  to  compute 
CCC  and  between  Q  and  AT,  and  the  coupling  coeffi¬ 
cient  for  Q  =  kAT,  at  day  d,  which  is  computed  in  this 
paper  from  a  91 -day  subset  between  45  days  prior  to 
and  45  days  after  that  date. 

The  global  distributions  of  (Rq^Ti)  for  the  30-day 
smoothed  data  (Fig.  8a),  the  60-day  smoothed  data  (Fig. 
8b),  and  the  90-day  smoothed  data  (Fig.  8c)  show  high 
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Gradient-type  without  solar  radiation 


month  /  date 


Gradient-type  without  solar  radiation 


month  /  date 


Gradient-type  without  solar  radiation 


month  /  date 


Fig.  9.  Seasonal  variation  of  coupling  coefficient  k  averaged  within 
five  geographic  zones,  determined  by  no  loss  of  variance  for  Q  = 
kAT2  after  the  daily  NCEP  dataset  undergoes  (a)  30-day  running 
mean,  (b)  60-day  running  mean,  and  (c)  90-day  running  mean.  Here, 
the  thin  solid  curve  indicates  the  southern  middle  and  high  latitude 
zone  (south  to  30°S),  the  thick  solid  curve  indicates  the  northern 
middle  and  high  latitude  zone  (north  to  30°N),  the  dotted  curve  in¬ 
dicates  the  southern  subtropical  zone  (10°-30°S),  the  dashed-dotted 
curve  indicates  the  northern  subtropical  zone  (10°-30°N),  and  the 
dashed  curve  indicates  the  equatorial  zone  (10°S-10°N). 


values  of  CCC  in  middle  and  high  latitudes  (higher  than 
30°)  away  from  coasts,  low  values  of  CCC  in  equatorial 
regions  (10°N-10°S),  and  not  much  difference  in  CCC 
among  the  three  smoothed  datasets.  Comparison  be¬ 
tween  Fig.  8  and  Fig.  lb  indicates  the  increase  of  CCC 
in  middle  and  high  latitudes  from  the  unsmoothed  data 
(~0.8)  to  the  smoothed  data  (—0.9).  Thus,  the  gradient- 
type  Haney  boundary  condition  is  also  valid  in  middle 
and  high  latitudes  (higher  than  30°)  away  from  coasts 
and  invalid  in  equatorial  regions  (10°N-10°S)  for 
smoothed  data.  The  validity  is  enhanced  when  the  data 
is  smoothed  over  30  days  or  longer. 

The  coupling  coefficient  k  between  Q  and  A Tl  for 
the  30-day  smoothed  data  (Fig.  9a),  the  60-day 
smoothed  data  (Fig.  9b),  and  the  90-day  smoothed  data 
(Fig.  9c)  has  large  seasonal  variations  in  the  equatorial 
and  subtropical  zones,  and  small  seasonal  variations  in 
the  middle  and  high  latitude  zones.  This  feature  is  very 
similar  to  the  unsmoothed  data.  Comparison  between 
Fig.  9  and  Fig.  6b  shows  no  drastic  change  of  k  from 
the  unsmoothed  daily  dataset  to  three  different  (30-day, 
60-day,  and  90-day)  smoothed  datasets  for  both  northern 
and  southern  middle  and  high  latitude  zones.  We  may 
use  70  W  m  2  K  1  for  the  northern  middle  and  high 
latitude  zone  and  65  W  m  2  K  1  for  the  southern  middle 
and  high  latitude  zone. 

8.  Conclusions 

Haney-type  surface  thermal  boundary  conditions  con¬ 
nect  net  downward  surface  heat  flux  to  air-sea  temper¬ 
ature  difference  (gradient-type  condition)  or  to  climate/ 
synoptic  sea  temperature  difference  (restoring-type  con¬ 
dition).  On  the  basis  of  cross-correlation  and  variance 
analyses  on  the  NCEP  net  downward  surface  heat  flux 
and  air-sea  temperature  data  during  1  October  1994- 
31  December  1995,  we  obtained  the  following  results: 

1)  The  restoring-type  conditions  do  not  represent  the 
surface  thermal  forcing  anywhere  in  the  world 
oceans. 

2)  For  the  equatorial  and  subtropical  oceans,  the  gra¬ 
dient-type  conditions  are  not  good  approximations 
for  the  surface  thermal  forcing. 

3)  For  the  middle  and  high  latitudes  away  from  coasts, 
the  gradient-type  conditions  are  good  approximation 
for  the  surface  thermal  forcing.  This  is  based  on  the 
high  correlation  between  net  downward  heat  flux  and 
air-sea  temperature  difference  and  associating  quasi¬ 
steadiness  of  the  coupling  coefficient  k.  Further¬ 
more,  there  is  a  better  correlation  when  the  solar 
shortwave  component  is  treated  separately. 

4)  A  value  of  70  W  m-2  K  1  (65  W  m-2  K1)  for  the 
coupling  coefficient  k  is  suggested  for  northern 
(southern)  middle  and  high  latitude  zones,  no  matter 
whether  the  data  is  smoothed  or  unsmoothed.  The 
suggested  values  are  about  twice  those  generally 
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used  (10-50  W  m  2  K  1 ).  This  might  increase  the 
net  air-sea  heat  flux  and  shorten  the  relaxation  time. 

We  also  tried  a  coupled  air-ocean  model  developed 
by  Russell  et  al.  (1995),  and  we  found  essentially  similar 
results. 
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